Collective shuttling of attracting particles in asymmetric narrow channels 
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The rectification of a single file of attracting particles subjected to a low frequency ac drive is 
proposed as a working mechanism for particle shuttling in an asymmetric narrow channel. Increasing 
the particle attraction results in the file condensing, as signalled by the dramatic enhancement of 
the net particle current. Magnitude and direction of the current become extremely sensitive to the 
actual size of the condensate, which can then be made to shuttle between two docking stations, 
transporting particles in one direction, with an efficiency much larger than conventional diffusive 
models predict. 

PACS numbers: 05.40.-a, 05.60.-k, 82.70.-y 



o 

(D 

e 

I 



o 



> 

cn 
in 

o 
o 



X 



Diffusion and directed transport of finite size interact- 
ing Brownian particles in constrained geometries are is- 
sues of growing relevance in soft matter physics In 
particular, the dynamics of one-dimensional (Id) gases, 
or single files, has found wide applications, ranging from 
the modeling of natural processes, such as in molecular 
motors and ion channels [2], to the design of artificial 
devices, involving, e.g., engineered porous materials Q, 
colloids [1, Q , and superconducting vortices ^ . 

A key feature, which is intrinsic to Id files in both 
natural and artificial channels, is attractive particle in- 
teractions, whose origin has been variously tracked down 
to excluded volume effects [2], Van der Waals forces 8], 
or the presence of additional passive molecules, like in 
colloid-polymer mixtures [§] . It has been recognized that 
pair attractive forces can induce the formation of large 
molecular clusters and that a dramatic increase in the 
overall file diffusion and mobility can ensue as a non- 
commensuration effect of condensate size and channel 
spatial periodicity [lo| . More generally, the incommen- 
surability between the lateral dimensions of a molecule 
and a catalyst is responsible for the "shape selectivity" 
of many a biological channel [llj . 

In this Letter we combine three different and important 
mechanisms, namely, (i) condensation in single files, (ii) 
Id transport and diffusion, and (iii) the ratchet effect and 
its applications, to outline the blueprint for a new class 
of nanoparticle shuttles. We show that Brownian parti- 
cles moving along a narrow channel with a ratchet-like 
profile, subjected to a low frequency center-symmetric 
ac drive, can condense to a cluster when the tempera- 
ture is lower than a certain threshold. This results in 
the net time-averaged particle current (J) being much 
larger (up to over one order of magnitude) than in the 
presence of weak or no attraction, when the particles 
are evenly distributed throughout the channel. More- 
over, (J) strongly depends on the ratio hN/L, where h 
is the particle length, N is the number of particles in 
the system, and L is the ratchet spatial period. Thus, in 
contrast to earlier reports, the ratchet current through 



the channel not only is strongly enhanced by condensa- 
tion, but it can also be inverted by adding or subtracting 
particles to the cluster (for appropriate h/L, just one!). 
This suggests a new concept of collective particle shuttle, 
where the motion of the cluster gets reversed any time, 
say, one particle is unloaded (uploaded) at a docking sta- 
tion (respectively, a sink or a source). Besides applying 
to the diffusion of long molecular chains in zeolites and 
the shape selective control of catalytic reactions in living 
cells, collective shuttling can also answer the longstand- 
ing question as why proton and electron biological pumps 
work much more effectively than predicted by purely dif- 
fusive shuttle models p^ . 

To demonstrate how ratchet effect can be controlled 
in the presence of attracting particles, we study a model 
system consisting of N Brownian particles moving along 
a circular Id channel. The particles interact via the pair 
potential 'w{xij) = w-hr{xij) + Wat{xij), which consists of 
a hard-core repulsive, Whr, and a long range attractive 
part, Wat, with Xij — \xi — Xj\ denoting the distance 
between the centers of particles i and j. In the following 
we assume Wi,t{xij) = — aexp(— Ax^j), where A^^ and 
a characterize, respectively, the range and the strength 
of the pair attraction. The dynamics of the particles 
is governed by a set of coupled overdamped stochastic 
Langevin equations (in rescaled units) 
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Here (^iit) is a stochastic white noise with zero mean 
and autocorrelation function {^i{t)^i{t')) — 6{t — t'), 
T is the temperature, and ^{{xj},t) is the total po- 
tential energy of the system, $({a;j},t) — ^i[U{xi) — 
F{t)xi] + \ J2j^i J2i wi\xi - Xj\), where U{x) represents 
the ratchet potential and F{t) an external ac drive. We 
choose a biharmonic function with period L = 1 for the 
ratchet potential, U{x) = sin (27rx) + 0.25 sin (47rx) [see 
Fig. mb)], and a square- wave with amplitude A for F{t). 
Finally, we set periodic boundary conditions at the end- 
points of the channel, the length of which is a multiple 
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oi L, S = ML. In doing so, we implicitly assume that 
hN < ML and that Wat{x) effectively vanishes for dis- 
tances comparable with S, w{S) ~ 0. 

To investigate the non-equilibrium properties of the 
system we numerically integrate the Langevin equations 
([!]). For computational reasons, in our simulations we re- 
place the hard-core potential w^r by the soft-core power- 
law potential Wsixij) = e{h^/xij)^^. For fixed e and /i* 
the effective hard core length hcs of these particles is a 
function of a and T (isj : in our simulations e = 0.01 
and h^, — 0.2 correspond to hgg ~ 0.16 for a = 10. We 
also study the system in the framework of the dynamical 
density functional theory (DDFT) pj]. DDFT approx- 
imates the Fokker-Planck equation for the one-particle 
density distribution p{x, t) to: 



dp{x,t) 
dt 



d_ 

dx 



p{x,t) 



d SF[p{x,t)] 
dx Sp{x,t) 



(2) 



where F[p\ is the Helmholtz free energy functional ob- 
tained from the equilibrium density functional theory 

[T6j . For hard rods with no attraction, Wat — 0, 
an exact analytical derivation yields the equilibrium 
free energy term i^hr = T j°°^dx p{x,t){\n p{x,t) — 
1) + UZodx^[p{x,t)Mx + h/2,t) + p{x - h/2,t)], 
with 0[/5(a;,i)] = — Tin [1 — 77(2;, i)] and ri{x,t) ~ 
J^^l^l^ dy p{y,t). Next, by adding a mean-field term 
to account for the particle attraction [l^, we obtain 
F[p\ = Fhr + fat, where F^t = ^ J dx J dywat{\ x - 
y \)p{x)p{y). Moreover, in the present case the func- 
tion p(a;, t) in Eq. ([2]) obeys periodic boundary condi- 
tions in space with period equal to the system size, 
namely, upon setting the center of the system at x = 0, 
p{-S/2,t) = p{S/2,t). Finally, on rewriting Eq. © as 
dp{x,t)/dt = —dJ{x,t)/dx, we introduce the instanta- 
neous particle current density, J{x,t). In the presence 
of the dc drives F = ±A, the unidirectional currents 

are computed as the time average of the instanta- 
neous current divided by the number of particles N, i.e. 

J± = lmit^^{l/Nt) dr J%^^^J{x,T)dx. To ex- 
tract from Eq. ([2]), we perform a Fourier mode ex- 
pansion of p{x, t) and then apply a standard continuation 
technique jl8l |. 

Two transport modes. When the pair attraction is suf- 
ficiently strong, the particles in the channel condense to 
form a compact cluster, which behaves like a single com- 
posite object. This phenomenon may be understood by 
considering a free single-file with F{t) = and U{x) = 0. 
The equilibrium properties of such a homogeneous sys- 
tem are well known in the thermodynamic limit 
{N,S — > 00 with constant density p = N/S). As the 
attraction strength, a, is increased, the particles tend to 
group together. Above some threshold value of a, in 2d 
and 3d systems the homogeneous state undergoes a phase 
transition to a state characterized by the coexistence of 
a low-density gas and a high-density liquid phase. In Id 
systems, particles also tend to cluster together on raising 
a, but no true phase transition can be defined [loj . 
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FIG. 1: (Color online) (a) Phase diagram of a infinitely long 
single-file with A = 3, T = 1, /i = 0.2. In the shaded region 
configurations with uniform density are linearly unstable due 
to particle attraction, (b) Density profiles p^{x,t) in a finite 
system of length S* = 10 and with N — 5 particles, obtained 
by solving Eq. ((2}. The dashed line represents the diffusion 
dominated case at a = 0. The density profiles labeled ti-t4 
are four snapshots of the traveling condensed mode (the arrow 
indicates the direction of motion) for A — 1, a — 10 and the 
remaining parameters as in (a). U{x) is also displayed for 
convenience, (c) Instantaneous current J{t) (solid line) and 
form factor M{t) (dashed line) vs. time for the simulation 
parameters of panel (b) . The vertical dashed lines denote the 
snapshot times ti-t4 also used in (b). 



Within our mean-field DDFT treatment, the homoge- 
neous one-body density p{x) becomes linearly unstable 
when the free energy per volume F[p]/S turns concave, 
i.e. for S^F/Sp^ < 0. As a consequence, inside the shaded 
(spinodal) region of Fig.[TJa), the attractive interactions 
are strong enough to overcome the stabilizing action of 
diffusion, thus making an initially homogeneous single- 
file unstable. For a system with finite size S and finite 
number of particles TV, such an instability signals the 
transition from an homogeneous state with evenly dis- 
tributed particles, to a condensed state characterized by 
a single cluster of approximate length hN and diffusing 
with effective temperature T/N. 

In the presence of periodic pinning U{x) and non-zero 
dc drive, condensation explains the emergence of two dif- 
ferent transport modes: a stationary homogeneous mode 
and a condensed traveling mode. To illustrate the tran- 
sition between these two modes, we set T = 1, A = 1 
and gradually increase a. For small a, the density pro- 
files p^{x), corresponding to F — ±1, are pinned by 
U{x), namely, are stationary functions with spatial pe- 
riod L = 1. This means that the particles are evenly 
distributed among the local minima of U{x), jumping 
occasionally from one well to the next in the direction of 
the drive, as illustrated in Fig.[Trb) for F = \. This is 
the standard ratchet mechanism [l| , which dominates at 
high temperatures or, equivalently, for weak interactions 



3 




FIG. 2; (Color online) (a) Unidirectional currents J vs. a, 
for T = 0.5, A = 1, iV = 5, Af = 5, A = 3, and ft = 0.16. The 
solid and dashed lines are, respectively, the stable and un- 
stable solutions of Eq. ((2]); symbols represent the correspond- 
ing simulation data. The Hopf bifurcations are marked by a 
square labeled HB. J* for h = are displayed for compari- 
son (dotted lines), (b) Zero-frequency rectification currents, 
(J), obtained from the data in (a). The positive a — )■ oo 
asymptotes coincide with the point hN — 0.85 in Fig. [Sjc). 



or dilute single-files. On increasing a this mode becomes 
unstable and is eventually replaced by a condensed trav- 
eling mode, resembling a wave packet of width hN trav- 
eling in the direction of the drive. Different phases of 
the traveling mode are shown in Fig.[TJb) at four differ- 
ent time instants, ti < t2 < < t^. As the N = 5 
particle cluster moves over the barrier of the pinning po- 
tential U{x), the interparticle distance gets modulated in 
time, giving rise to spatio-temporal cluster oscillations. 
The period of these oscillations coincides with the time 
needed for the cluster to move across a unit cell of U{x). 
This barrier crossing mechanism is better visualized in 
Fig.[IJc), where the instantaneous velocity of a cluster 

particle, J{t) ^ {^/N) J_g/2 Ji^^t) d-^^ is contrasted with 
the form factor, M (t) = ptiax~ Pmin^ introduced to quan- 
tify the spatial oscillations of the cluster. Here p^ax i^ 
the absolute maximum and p^j^^ the largest of the local 
minima of the instantaneous cluster density p^{x^t). 

Current enhancement. The efficiency of particle trans- 
port along the channel is characterized by the unidirec- 
tional currents driven by -F = ±A. Correspondingly, 
the rectification current, (J), induced by a square- wave 
ac drive F{t) of amplitude A and vanishingly low fre- 
quency, is simply the average (J) = \{J^ + J^)- In 
Fig.UJa) we compare our DDFT results for the two cur- 
rents J* as a function of a (solid lines) , with the simu- 
lation outcome (symbols). As expected, for small a the 
long time limit solutions of Eq. ^ are stationary func- 
tions with spatial period L. The branch of the L-periodic 
stationary solutions becomes unstable via a Hopf bifur- 
cation (HB) at auB — 2; see Fig.[2Ua). For a above the 
HB threshold auB, a new branch of stable time-periodic 
solutions appears. The transition from the stationary to 
the time-periodic regime corresponds to the transition 
from the stationary to the traveling condensed mode. At 
the onset of the traveling mode, immediately to the right 
of the HB point, both | increase sharply with a, well 
above their zero attraction values, until, for a « 5, they 



approach the expected upper bound, | = A (unpinned 
single- file regime [20l|). 

The transport enhancement in the traveling mode is 
not solely a consequence of the particle attraction: by 
preventing condensation to a vanishingly small volume, 
hard-core particle repulsion also plays a crucial role. If 
one sets = 0, no enhancement of the particle mobility 
with increasing a occurs [dotted lines in Fig.l^^a)]; on 
the contrary, for pointlike particles drop by as much 
as one order of magnitude after condensation sets in, as 
already pointed out in Ref. [2]| . 

The difference between the unidirectional currents 
depends both on the interaction strength a and on the 
particle length h. Indeed, the a dependence of the zero- 
frequency ratchet current, (J) = ^(J+-|-J~), may widely 
vary with h. In the example of Fig.[5Ub), |(J)| increases 
rapidly, retaining its initial [a — 0) negative sign. How- 
ever, for a > 5 it reverses sign and eventually levels off 
with even larger a modulus. Note that the existence of a 
zero-crossing in the curve (J) versus a is determined by 
the sign of the curve asymptote for a — >■ cx). In this limit 
(J) oscillates with h as one better sees by investigating 
our model in the strong attraction approximation. 

Cluster size and current reversals. In the limit a — )■ oo, 
when particle attraction dominates over thermal diffusion 
and static pinning, Eqs. ([T]) can be reduced to a single 
equation describing the evolution of the center of mass, 
y — {l/N) Xi, of the particle condensate, 

dy _ Ujy + hN) ~ U{v) 
dt hN 

where ^(i) has the same statistics as ^i{t) in Eqs. ([T]). 
This equation helps clarify the role of the particle size, 
h, and the number, N, of particle in the cluster in the 
depinning of the condensed mode. By inspecting Eq. ([3]), 
we notice that for large N diffusion becomes negligible, 
and the unidirectional currents set in only when the 
drive amplitude A overcomes the pinning force exerted 
by U{x), either from the right, Fr, or from the left, Fl 
(both Fb,,l are functions of hN). Note that for a sin- 
gle particle, N = 1, or for hN 0, F^i — 3tt and 
Fl — 37r/2 [Ij. Upon neglecting thermal fluctuations, 
and the zero-frequency rectification current, (J), can 
both be computed analytically from Eq. ([3]) [l| as func- 
tions of hN. For hN and F chosen within the shaded re- 
gion of Fig.[3l^a) the condensed mode is pinned. Nonzero 
currents thus emerge only due to thermal fluctuations 
with effective temperature T/N. Most remarkably, for 
appropriate combinations of h and N, one can achieve 
complete locking of the condensed mode in one direc- 
tion, but not in the other [2l| , leading to the observation 
that for a certain value of A the modulus of (J) cannot 
exceed Jmax = A/2. This is a benchmark against which 
to compare the efficiency of single-file rectification. 

For finite T, depend on h and N separately. In 
Fig.[2Ib) we plot versus N in three different regimes: 
in the a ^ oo limit described by Eq. ([3]) with h = 0.16 
(solid lines); from numerical simulation of the original 
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FIG. 3: (Color online) (a) Pinning {F, hN) region (shaded) 
for the condensed mode in the limit a — >■ oo, Eq. ([3]), and 
T = 0. The horizontal lines correspond to _F = ±1 used in 
the simulations, (b) Solid lines: J* vs. hN, obtained inte- 
grating Eq. ((31 for h = 0.16; symbols: data points obtained 
via simulation of Eqs. ^ with a — 10, A = 3, hes — 0.16 and 
varying A'^; dashed lines: J* vs. hN computed from Eq. ([2]) 
with a — 0. Other parameters are A = 1, T = 0.5, and 
M = 10. (c) The zero-frequency rectification currents (J) 
vs. hN computed using the data in (b). A qualitatively sim- 
ilar plot is obtained on varying h for fixed A^. The current 
peaks should be compared with Jmax ~ 0.5. (d) Sketch of the 
molecular shuttle described in the text. 

system, Eq. ([T|) with strong pair attraction, a — 10, 
A = 3, M = 10 and hcs ~ 0.16 (symbols); and from the 
DDFT approximation, Eq. for a = 0, h = 0.16 and 
M — 4 (dashed lines). All other parameters, T — 0.5, 
A = 1, are the same. The corresponding zero-frequency 
rectification currents, (J), are displayed in Fig. |3jc) . Sim- 
ilar behaviors can be obtained by varying h at constant 
N. In conclusion, on increasing the size of the condensed 
cluster our model exhibits: (i) an oscillatory growth of 
J^, which eventually tend to the respective (unpinned 



particle) limits = ±1, and (ii) sign reversals of (J) 
for hN an integer multiple of L/2. The qualitative ex- 
planation of this behavior is simple for a — > oo, when the 
single- file of particles can be regarded as a continuous rod 
of length hN. The total pinning force exerted by U (x) on 
such a traveling rod vanishes for hN — nL, with n a posi- 
tive integer (zero effective pinning) , whereas its intensity 
is maximum for hN w (2n — l)L/2. Correspondingly, 
= ±1 for hN an integer multiple of L [Fig.[3fb)] and 
(J) = for hN an integer multiple of L/2 [Fig.[3Kc)]. For 
finite a the effective rod is somewhat larger than hN, so 
that zeros of (J) may occur for hN < 1, like in Fig.[5l^b). 

Collective shuttles. Particle size and number can be 
used to selectively control the sign and magnitude of the 
net current. This concept is illustrated and applied in 
the toy molecular shuttle, sketched in Fig. Eld), where 
a cluster of = 2 rod-shaped molecular units, each of 
length h = 0.4 (in dimensionless units), moves to the 
right [Fig.[3[c)] until it docks at the loading station at 
the exit of the channel (source) . After binding to an ad- 
ditional molecular unit, the enlarged cluster now reverses 
its drift velocity [Fig.|3Uc)], thus dragging its cargo to the 
unloading station at the left exit of the channel (sink), 
where it releases the leftmost unit, so that the process 
can repeat itself. A steady molecular flux is thus estab- 
lished along the channel between the source and the sink, 
more efficiently than by mere single molecule (or, worse, 
single-file) diffusion. Consequently, our shuttle concept 
provides a much more efhcient mechanism for ion trans- 
port in bio pumps than standard diffusion shuttles jl^ . 
Finally, we stress that high-performance collective shut- 
tles only work with entrained particles of finite-size; for 
pointlike particles the resulting condensate cluster would 
have zero width and the file current would be suppressed. 
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